Sequential Implementation of Global Quantum Operations 
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We study the possibility for a global unitary applied on an arbitrary number of qubits to be 
decomposed in a sequential unitary procedure, where an ancillary system is allowed to interact only 
once with each qubit. We prove that sequential unitary decompositions are in general impossible for 
genuine entangling operations, even with an infinite-dimensional ancilla, being the controUed-NOT 
gate a paradigmatic example. Nevertheless, we find particular nontrivial operations in quantum 
information that can be performed in a sequential unitary manner, as is the case of quantum error 
correction and quantum cloning. 
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Unitary evolution, a fundamental pillar of quantum 
dynamics, can be described as the action of a time- 
dependent unitary operator on an initial quantum state. 
In quantum information, the resulting evolution of a 
global quantum operation acting on N qubits can be 
rewritten as a quantum circuit made out of one- and 
two-qubit gates. The application of these gates is nonse- 
quential in general, and the same qubits can be operated 
repeatedly at different stages of the prescribed quantum 
algorithm ^l| . An arbitrary global unitary is able to pro- 
duce a unitary transformation that maps N qubits to N 
qubits. In cases where M of the initial N qubits are un- 
known and the rest are blanks, as happens in quantum 
cloning Q and quantum error correction protocols 
an M — » iV qubit map that is formally described by an 
isometric operation can be considered. An isometry is a 
unitary operation acting only on a subset of its domain. 

Engineering a realistic arbitrary global unitary U act- 
ing simultaneously on N qubits is known to be a diffi- 
cult problem. It may be then desirable to decompose U 
into a sequential unitary procedure, where each qubit is 
allowed to interact locally and only once with an itiner- 
ant ancillary system A, and without measurements (see 
Fig. 1). The picture of a sequential quantum factory 
where a global U is decomposed into a sequence of lo- 
cal qubit-ancilla unitary steps is a valid one and could 
be related to particular quantum Turing machines 
However, one should be careful to differentiate this dis- 
tinct physical problem with the already studied sequen- 
tial quantum factory of states [1, [6|. The latter was 
shown to be always possible for any state, as long as 
the dimension of the ancilla matches the bond dimension 
of its corresponding Matrix Product State (MPS) repre- 
sentation [tI . It is noteworthy to mention the unexpected 
connection to MPS theory, which has recently shown its 
power for understanding spin chain correlations [Sl, clas- 
sical simulation of quantum entangled systems [9|], and 




FIG. 1: Sketch of the proposed sequential quantum factory 
of operations. When possible, it decomposes a global uni- 
tary operation U acting on A'^ qubits into a sequence of lo- 
cal unitary steps where each red qubit interacts only once 
with the ancilla A, and measurements are not allowed. When 
each qubit-ancilla interaction is finished, the corresponding 
red qubit in the chain turns into green and the following red 
qubit starts its interaction with A. 



density- matrix renormalization group techniques [10| . 

In general, sequentiality and the absence of measure- 
ments may reduce enormously the complexity in the 
physical implementation of a global unitary, though they 
are not unavoidable requirements. However, different ex- 
periments and physical settings demand or are restricted 
to sequential operations for achieving specific quantum 
information tasks, and allowing measurements may lower 
the quality of deterministic protocols. Paradigmatic 
examples are the sequential generation of multi-qubit 
photonic chains, entangled in energy or polarization 
states [llj , and the sequential generation of electron- 
spin entangled states in quantum dot setups (l3 |. In 
the former example, the ancillary system can be a cav- 
ity mode, while in the latter it can be a mixed nuclear 
spin bath. We point out that our proposal is related to 
the quantum interface problem [ll|, studying the faith- 
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ful transfer of information from one quantum node to 
another. From the point of view of unitary operations, 
Delgado et al. proposed recently a sequential imple- 
mentation of quantum cloning from 1 to A'^ copies, but in 
their scheme qubit measurements were allowed to achieve 
the desired goal. The case of one-way quantum comput- 
ing 1J| is special, it is deterministic and sequential, but 
it is measurement based and requires the generation of 
a particular initial entangled state. Below, when we re- 
fer to a sequential decomposition we assume that it is a 
unitary one, unless other specification is given. 

Strikingly, we will prove below that a sequential de- 
composition of a global unitary U acting simultane- 
ously on N qubits is impossible in general, even with an 
infinite-dimensional ancilla. By using MPS theory, we 
will also show how to discriminate the global operations 
that can be decomposed in a sequential manner and, in 
the affirmative case, provide with the optimal dimension 
of the ancilla. We will illustrate our findings with the 
Shor nine qubit code [15,] and the optimal 1 ^ N cloning 
machine [21 , showing that the optimal dimension needed 
for the ancilla in both cases is rather small. 

We first illustrate this surprising result with such a 
simple example as the paradigmatic two-qubit controUed- 
NOT (CNOT) gate, which flips the target qubit whenever 
the control qubit is in state |1) and leaves it unchanged 
when the control qubit is in state |0). We consider a 
CNOT gate, where qubits 1 and 2 are the control and 
target qubits, respectively, acting on the initial tripartite 
state I a) 1 1 6) 2 1 c) a. The latter element in this tripartite 
product state represents an ancillary system A with ar- 
bitrary dimension D. 

If we could decompose the CNOT sequentially, we 
could write 



CNOT12 <E>la = U^aUla- 

Therefore, the input | + 00)i2a, |+) being ^(|0} -I- 
would transform into --i=(|00) + |ll))i2|0)a. This implies 
(looking at the entanglement in the partition l|2a) that 
Uia\ + 0) has to be maximally entangled. However, | -t- 
+0)i2o remains untouched by the CNOT, which implies 
(again looking at the partition l|2a) that Uia\ + 0) is a 
product state, giving the desired contradiction. 
We prove now the following generalization 
Theorem 1: No entangling unitary can be implemented 
in a sequential way. 

Let us call n the smallest i for which Uia- ■ - Uia = 
Ui...i ^ C/o, that is, the smallest step in which the ancilla 
always decouples from the qubit chain. This number ex- 
ists since it happens at least in the A^-th step. To finish 
the argument it is enough to show that n = 1. If it were 
not true, we could group the first n — \ spins and the n-th 
one and then have U2aUia = U12 ® la, where Uia entan- 
gles 1 and a (by the choice of n). Therefore, assuming 
without loss of generality that J7ia|00) is entangled, we 
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FIG. 2: Trivial protocol for the sequential decomposition of a 
1 A'^ isometry. First, the desired isometry is implemented 
locally between qubit 1 and A'' — 1 qubits in the ancilla degrees 
of freedom. Second, each ancilla qubit is mapped locally to 
their corresponding blank chain qubits. 



may write 



0= (la|C/2af/la|0j0)i2a 



El 

4=0 



{la\U2a\^J)a2, (1) 



with r > 1. Since the are linearly independent, 

{'i-a\U2a\ij)a2 = for all j and i = 0, 1, implying that 
det([/2a) = and producing the desired contradiction. 

It is noteworthy to mention that, in contrast to the pre- 
vious negative result, one can always implement sequen- 
tially a 1 ^ iV isometry as long as the ancilla dimension 
is large enough, see Fig. 2. The way to see it is straight- 
forward: i) we apply locally the desired global isometry 
on qubit 1 and N — I blank ancilla qubits, ii) we imple- 
ment local swap gates between the ancilla qubits and the 
A'^ — 1 blank chain qubits. Clearly, the worst case scenario 
will require an exponentially large ancilla dimension but 
this fact does not invalidate the previous general claim 

Therefore, the remaining task is to find the optimal 
dimension of the ancilla needed to implement a given 
1 ~> N isometry. The solution to this problem will come 
via the use of MPS 0] representations of vectors. They 
can be defined as 



E 



IN]^[N-1] 



■AfJAf^\^l^2■■■^N). (2) 



where A ™' are Dm+i x Dm matrices with Di = D^+i ~ 
1. We define the MPS bond dimension by D = max^ Di 
and it can be shown 0, 0] that there exists a canonical 
form for any MPS that gives the minimal bond dimen- 
sion. This MPS canonical form is defined by the condi- 

tions: (i) E.^"''^!"' = 1, (ii) 4"' M^-iUl™!^ = 
A['"l, for all 1 < m < iV, and (ui) AI"! = Al^l = 1 with 
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FIG. 3: Diagram showing the sequential decomposition of a 
1 ^ A'^ isometry U, acting on an arbitrary initial state \^p}. 



each A[™1 diagonal, positive, full rank and tiAl™! = 1. 
Since it is obtained by performing successive Schmidt de- 
compositions, this canonical form is essentially unique. 
Moreover, using Theorem 3.1.1' of [l^l, one can see that 
all canonical forms are related by the condition 



(3) 



with Vb = Vjv = 1 and Vm unitary such that SVm , Am] = 
for all m. 

In our case we want to get a canonical MPS represen- 
tation of the 1 ^ iV isometry IJ . We do it by artificially 
grouping the input and output of the first qubit to get 



U = 



E 



Jl,Jl,...,ljV 



(4) 



The key observation is to note that the relations 

ii ,k 

allow us to identify the existence of a sequential decom- 
position of U, as in Fig. 3, with an MPS decomposition 
of U as in ([4]) satisfying the additional properties 



= 1, m > 2. 



(6) 

Our main result states that the optimal ancilla to im- 
plement U sequentially is exactly the hand dimension of 
its canonical MPS form. Moreover, any of these canon- 
ical MPS representations gives via 0^ the sequential de- 
composition of U as in Fig. 3. 

The proof of this result is based on recent mathemat- 
ical studies in MPS theory 0, Q . We start then with a 
sequential decomposition of our given 1 — > isometry, 
for instance the one given by Fig. 2. Using ^ this gives 



an MPS representation for U , using matrices B. We will 
show that from there one can obtain a canonical form for 
U , with matrices A, and still satisfying Using ([3]) 
this in turn implies that all canonical forms for U satisfy 

The procedure to go from matrices B to A was de- 
veloped in the proof of Theorem 2 in Perez-Garci'a et 
al. [7|. There it is shown how to construct in general 
non-square matrices Zj with Z^^Zj = 1 such that, if we 



define A™.. = ZJsW^-^ , ^[^1 



Bfjz^^, and 4 



■ii ji J- n ji ' «iv 
Zl^B^"^ Z„i-i for 1 < m < A^, the A's give a canonical 

form for the U operator. Moreover ZiZli?]^',, = s!^' ■ . 
But we finish now by noticing that 



A™., = y Sflt z.zlB^'K, 



(7) 



(8) 



A useful consequence of the previous result is the fol- 
lowing 

Corollary: If Di is the optimal bond dimension for 
J7|i), then the optimal ancilla dimension to generate V 
sequentially is upper bounded by Do + D\. 

To see it, let us take MPS representations of ?7|0) and 
?7|1), given by matrices A and B respectively. Since 
V = ?7|0)(0|-|-C/|1}(1|, to finish the argument it is enough 
to notice that one has the MPS decomposition U — 



[N] 



c 



•4^5J*i)0'iI ® h---^N), taking 
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B 



[1] 

41,0 







B. 



[1] 



(9) 



(10) 



Let us illustrate the strength of the result with a cou- 
ple of important examples. The first one is the nine- 
qubit code of Shor [l5[, which was the first quantum 
code that allows to correct an arbitrary one-qubit er- 
ror. It encodes one qubit into nine, via the encoding 
map |0} ^ (|000) + |111))®3, |1) ^ {\000) - |111))®3. 
Since, using the notation of the Corollary, Dq = Di = 2, 
one can implement this encoding operation sequentially 
with an ancilla consisting of only 2 qubits (dimension 4). 
This makes the sequential implementation of quantum 
operations appear as an interesting tool for fault toler- 
ant quantum computation. 

The second example is related to optimal quantum 
cloning In this context [17], one tries to overcome 
the restrictions imposed by the no-cloning theorem llSlI 
at the price of lowering the cloning fidelities. In Ref. [l9|. 
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it is shown how to clone in an optimal way, that is, 
how to construct a 1 ^ (27V — 1) isometry U such that 
for any of the first N qubits, the clones, one has (i) 
Pi = . . . = Pn = p where pi is the reduced density matrix 
for qubit i of U\il}) (all clones are equal), (ii) the fidelity 
{tp\p\'\j]) between the clones and the original state is max- 
imal. For this U it is shown in [l^ how Do = Di — N . 
Then, using the Corollary, we conclude that the optimal 
cloning isometry U of one qubit among N participants 
can be done sequentially with an ancilla of dimension 2N, 
that is, consisting of log(iV) + 1 qubits. 

We have already demonstrated that in the case of 
a full N —>■ N unitary, no entangling global unitary 
is decomposable in sequential qubit-ancilla local steps, 
whereas every 1 ^ isometry has a sequential decom- 
position. What happens then in the intermediate case of 
an M ^ iV isometry (M < N)? As one could expect, 
there are some nontrivial (entangling) unitaries that do 
admit a sequential decomposition. To fully characterize 
which is the case, one can also apply here the MPS ar- 
gumentations used in the 1 iV qubit case to get the 
following result. 

An M ^ N isometry U can be implemented in a se- 
quential manner if and only ify^,, Al-'^'t / = <5,-, i' 1 
/or one ( and then for all by Eq. (0) ) canonical form and 
forl<k< M. 

As a consequence, if U can be implemented in a se- 
quential manner, the optimal of such decompositions, in 
terms of the dimension of the ancilla, is given by any of 
the MPS canonical forms for U, written as 



U 



Ji, ■■■,JM ,^l,■■■,^N 



•••4'!,j*i---*^)o-i---jm|. (11) 

The relations between the sequential decomposition of 
U as in Fig. 3 and the MPS form (fTTj) are given now by 

^W|0)a|jl)=E(^^l4L|fc)aNl), 

ii,k 

V['=l|r.)abfc) = E i^k\A^L\^k)\s,)a\^k) , k<M 

^'"kfc)a|0) = Y.^-"^\^tk\^k)\sk)a\lk) , k>M. 
ik .Sk 

In conclusion, we have analyzed the sequential de- 
composition of arbitrary unitary operations and we have 
shown that this is unattainable in general for the N ^ N 
qubit case. On the other hand, we have proved that 
the 1 ^ A^ qubit case is always possible and we have 
found the optimal ancilla dimension. For the interme- 
diate M ^ N qubit case, with AI < N, we have given 
necessary and sufficient conditions for the sequential de- 
composition to be possible. We believe these results 
are of importance in the context of quantum theory and 



quantum information applications. With no doubt, they 
will prove to be relevant in any future implementation of 
quantum networks, quantum error correction, quantum 
cryptographic attacks, and related physical setups. 
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